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The critical point of a topological phase transition is described by a conformal field theory, where 
finite-size corrections to energy are uniquely related to its central charge. We investigate the finite- 
size scaling away from criticality and find a scaling function, which discriminates between phases with 
different topological indexes. This function appears to be universal for all five Altland-Zirnbauer 
symmetry classes with non-trivial topology in one spatial dimension. We obtain an analytic form 
of the scaling function and compare it with numerical results. 


Since the introduction of topological order in con¬ 
densed matter physics, the field of topological insulators 
received constantly growing attentiorPlEl. Although non¬ 
interacting topological phases were fully classifiec P^ a nd 
a plethora of topological edge states characterizecpEMII]^ 
little attention was given so far to finite-size effects 
around the topological transition. An important ques¬ 
tion is whether finite-size scaling is capable to distinguish 
between topological indexes and may be used as an in¬ 
dicator of the topological nature of the transition. One 
may also ask whether such scaling is universal or specific 
to a particular symmetry class, e.g. sensitive to Z vs. Z 2 
topological index. 

In this paper we discuss the finite-size scaling of the 
ground state energy across topological phase transitions 
in 1 -|- 1 dimensional models. The critical point in such 
models is described by a conformal field theorjff^ (CFT). 
The finite-size, N, scaling of the ground state e nergy 
E{N,0) for an open system at criticality is knowrP^*^^ 
to be 


E{N, 0) = N e-(0) + b{0) - ^ ^ (1) 

where e(0) is the average bulk energy per particle, 6(0) 
the size-independent boundary term and argument (0) 
specifies the exact critical point. Here length is measured 
in units of lattice spacing and energy in units of the Fermi 
velocity over the lattice spacing. The 1/A^ term appears 
to be universal and depends only on c - the central charge 
of the Virasoro algebrsP^. 

A relevant perturbation drives the system away from 
criticality, creating a spectral gap 2m and a correspond¬ 
ing correlation length ^ = 1/m. Our main observation is 
that the CFT expansion Q may be generalized as 

E{N, m) = N e(m) -|- 6(m) — ^ f{Nm) + , (2) 

where in the double scaling limifP^ A^ —>■ 00 and m —>■ 0, 
while w = Nm = N/^ = const, the function f{w), Fig.[^ 
is universal for all 5 Altland-Zirnbauer symmetry classes 
with non-trivial t opo logy in 1 spatial dimension (AIII, 
BDI, Dill, D, CIl)^. Hereafter we identify m > 0 with 
the topological and m < 0 with non-topological, or lesser 
topological index, side of the transition. Most notably. 



FIG. 1: (Color online) Numerical results for f{w), where 
w = Nm and N = 100 for 5 topologically non-trivial 
symmetry classes in one spatial dimension, (a) The case 
for open boundary conditions which is sensitive to 
topology, the solid line is the scaling function given by 
Eq. (b) In periodic boundary conditions the results 
are independent of the topological index and the scaling 
function Eq. (101 is symmetric. There is a difference 
between even (/(ui) negative) and odd (/(w) positive) 
number of unit cells N = 100,101. 
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the scaling function for open boundary conditions ex¬ 
hibits markedly distinct behavior on the two sides of the 
topological transition, while for periodic boundary con¬ 
ditions it is symmetricli^. Curiously, a similar scaling 
function for the entanglement entropjff^ appears to be 
symmetric across the topological transitioiP^, and it is 
only the iV-independent boundary term which is sensi¬ 
tive to the topological index. 

One may worry that in the double scaling limit depen¬ 
dence on m of the bulk and boundary terms should not 
be kept. This is not quite so, because of their singular 
dependence on the gap. As we explain below 

e(m) = e(0) — [m^ -I- C>(m^)] lna|m|; (3) 

b{m) = b{0) + — [m + 0{rri^)] lnQ;f,|m|, (4) 

where a and ai, are non-universal constants. As a re¬ 
sult the double scaling limit ([^ for the energy may be 
equivalently written as 

E{N,w) = Ne{0) + bi0)+^^^^{2w-w^)-^f2iw), 

( 5 ) 

where f 2 {w) = f{w) + logab|w| — log a|w|) in¬ 

corporates non-universal terms ^ w and ^ uP'. Since 
these latter may be easily subtracted both numerically 
and analytically, it is preferable to use the expansion Q 
with the fully universal function f{w). 

Universality of the scaling function: Before discussing 
analytic properties of the scaling function f{w) let us fo¬ 
cus on our numerical setup and demonstrate the univer¬ 
sal behavior for different symmetry classes. We consider 
models for all five Altland-Zirnbauer symmetry classes 
which are topologically non-trivial in one dimensiorP. To 
extract the scaling function we use Eq. ([^. E{N,m) is 
the sum of all eigenvalues which are obtained by numer¬ 
ical diagonalization. The average energy e is the integral 
over the entire Brillouin zone of the dispersion relation 
for all filled bands, which are calculated from the k-space 
representation of the Hamiltonian. The boundary term 
may be also calculated analytically (see below) or alter¬ 
natively approximated by b{m) « E{N,m) — Ne(rn) for 
some large N, say N = 1000. We have checked that the 
two ways are in excellent agreement. 

In Am symmetry class we use the standar d Su - 
Schrieffer-Heeger (SSH) tight-binding HamiltoniarP^l^oi 

N N-1 

'HaIII = '^tlc\jCB,j + ^ t2c\jCA,j+l + ^.C. ( 6 ) 

1=1 1=1 

Here we choose the gauge for the momentum to have the 
gap closing at fc = 0, where the dispersion relation reads 
£{k) = ±\/tl +t 2 — 2tit2 cosk. For ti p t 2 neighboring 
sites form dimers, where for ti > t 2 all sites are part of 
a dimer, but for ti < t 2 the two sites at the ends of the 
chain are unpaired (from now on m = t 2 — ti). Thus there 
are two distinct phases with topological index Z = 0 
OT Z = 1 respectively. In the case of n similar parallel 




FIG. 2: (Color online) SSH model (AIII symmetry class): 

(a) Visualization of the phase shift across the transition. 
For a fixed w = Nm, a state exists with energy P{k) if 
w = kN cot{kN), Eq. 0. In both limits w —>■ ±oo there 
are states at kN = mr, however one pair of states collides 
at ui = 1 (blue) and obtains imaginary k, shown as dashed 
line, (b) Energy spectrum near the gap for a A = 50 SSH 
Hamiltonian (black), with a pair of evolving edge states 
(blue). The edge state crosses the bulk Dirac cone (thin 
red) at UI = 1, where the momentum becomes imaginary, 
(c) Comparison of numerical results for SSH Hamiltonians 
with the scaling function 0 for system sizes 
N = 50,100, 200 and three transitions in a 2-chain SSH 
model with topological indexes Z = 0 —>■ 1,1 —>■ 2,0 —>■ 2. 


chains the topological index takes values Z = {0,1,..., n}. 
Fig# shows results for the scaling function at different 
system sizes, as well as three different transitions {Z = 
0 —1,0 —>■ 2,1 —>■ 2) in a system with two parallel 
chains. The results scale with w as the only parameter. 
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and all transitions agree with the analytic result, Eqs. ([^, 
The other symmetry classes are discussed in the 
supplemental material^. Here we only show numerical 
results in Fig.[2 which confirm universality across all five 
topological classes. 

Analytic properties: Universality of f{w) function is 
related to the fact that, similarly to the CFT result Q, 
it is fully determined by the vicinity of the critical point. 
One may thus approximate a near-critical system by the 
Dirac Hamiltonian, e.g. in AHI symmetry class, % = 
mai + idx<J 2 , where the Pauli matrices act in the space 
of A/B sublattices, cf. Eq. Q. Assuming that outside 
of the interval 0 < x < N the mass is very big and, e.g., 
negative one derives the boundary conditions 4'yi(0) = 
'I'b(A^) = 0. The quantized values of fc > 0 are given by 

TTi in 

cosikN5{k)) = 0] tanJffc) =-p =-r—. (7) 

k kN 

As a result the spectrum is determined by the condition 
w = Nm = kN cot{kN), plotted in Fig. and is given 
by e^(fc) = ±-\/m^ -I- k"^. At w = 1 two of its real so¬ 
lutions collide and switch to purely imaginary ones, cor¬ 
responding to the decaying edge states. Notice that the 
non-propagating states do not form at to = 0, as could 
be naively expected, but rather at to = 1/A^. 

Using the argument principle the total groundstate en¬ 
ergy is given by 

1 r 

E{N, m) = - (fc) dk In [cos{Nk + d(fc))], (8) 

for the dispersion relation e~{k) of the filled lower band. 
The contour runs in the complex A:-plane encircling all 
solutions of Eq. Q. The bulk and boundary terms are 
given by Ne + b = J{dk/2TT)e~ {k)[N + dk6{k)], where 
N + dkS{k) are bulk and boundary parts of the continu¬ 
ous density of states. To find the scaling function f{w) 
one subtracts Ne -I- b from Eq. (|^, deforms the integra¬ 
tion contour to run along the branch cut of Vto^ -I- 
and rescales the integration variable as z = ikN. As a 
result, one finds 


[ —\/z'^ — w'^dzln 

\ _|_ g-2z-25„(z) 

^\w\ 

L -1 


1 . 

where 6w{z) = — tanh (w/z). (For more detail on the 
derivation see the supplemental material^.) This expres¬ 
sion is plotted in Figs. as a full line and is in good 

agreement with the numerical data. 

Before discussing analytic properties of this scaling 
function let us add a couple of remarks: (i) though 
the derivation was given for the model in symmetry 
class AHI, the same logic works for the other symmetry 
classes. One needs to subtract proper model-dependent 
bulk and boundary parts, but the scaling term is only de¬ 
termined by the vicinity of the Dirac point and remains 
unchanged; (ii) a similar derivation may be applied to 
the case of periodic boundary conditions. In the gauge 


chosen after equation Q periodic boundary conditions 
give "^(N) = (—1)^5'(0). In this case the quantization 
condition 0 changes to cos{kN) = (—1)^. After sub¬ 
tracting the bulk energy (there is no boundary term in 
this case) and following the same steps one arrives at a 
similar scaling function: 


poo 1 

f{w) = —2 / — \/z'^ — w'^ dz In Tl — (— 

J\w\ 


l\w\ 

^ _2w ^ Kijjw) ^ -jv 

^ h J 


( 10 ) 


where Ki{x) is the modified Bessel function. This func¬ 
tion f{w) is manifestly symmetric across the topological 
phase transition, as it must be for periodic boundary con¬ 
ditions. However even in the scaling limit A^ —>■ oo it is 
dependent on parity of N, see also Fig.[^. The difference 
may be attributed to the level crossing at the gap closing 
point fc = 0 and w = 0 for even N, explaining ^ |w| non- 
analytic behavior of the scaling function. For odd N, all 
levels undergo avoided crossings and the scaling fun ction 
is free from such non-analyticity. The CFT resullP^*^] 
predicts /(O) = tt/6, which agrees with the case for odd 
N, while for even N we obtain /(O) = —tt/S. 

Returning to a system with open boundary conditions, 
at small |w| ^ 1 Eq. Q leads to: 

f{w) - ^ + ^{-2w + w^)\n\w\ + .... (11) 

The first term here is in agreement with the CFT limit 
Q. The subsequent terms ensure that f 2 iw) function, 
defined after Eq. ([^, is analytic. Indeed, at any finite N 
the ground state energy E{N, to) and all its derivatives 
must be non-singular at to = 0. To derive the second 
term in Eq. © one may employ monodromy transfor¬ 
mation, which rotates complex ic in a small circle around 
zercP^. Upon such transformation the right hand side of 
Eq. (§ picks up a contribution given by a closed con¬ 
tour integral around a branch cut — |rc| < z < |rc| times 
the number of revolutions. Calculation of such an in¬ 
tegral leads to i{—2w + w^), implying that f(w) must 
have logarithmic branch cut terminating at w = 0 with 
the discontinuity across it given by this value. Hence 
Eq. (111. We note in passing that in addition to such 


logarithmic branch cut, f(w) function has an inhnite se¬ 
quence of square root branch cuts along the imaginary 
axis of w. 

At large argument |w| 1, i.e. N |^|, the finite-size 

corrections decay exponentially. Remarkably the rate of 
the decay appears to be sensitive to the topology: 


(In Z symmetry classes, the two lines may be attributed 
to Z and ^ -|- 1 topological index.) The fact that on 
the topological side the scaling function decays half as 
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fast as on the trivial side may be associated with the ap¬ 
pearance of the edge states in the middle of the gap and 
effectively cutting the gap in half. In fact, the purely 
imaginary solution of w = kN cot{kN) at w 1 gives 
the energy of the edge states as e = ±2we~'"/N. This 
is identical to the asymptotic of f{w)/N on the topolog¬ 
ical side of the transition, Eq. (121, indicating that the 
latter originates solely from the edge state. In the case 
of periodic boundary cond itions the large w asymptotic 
is f{w) = — (—l)^-\/2/7r|r(;|^/^e“l’"l, which is different 
from both sides of the transition, Eq. (12), in the open 
boundary condition case. 

Eurthermore note that there develops a peak at w = 
iV/^ = 1 on the topological side (cf. Fig.[^. At this point 
there is a crossover between the regime of the correlation 
length being larger than the system size to smaller than 
the system size. In other words, here the two edge states 
at opposite ends transform from being delocalized and 
correlated to localized modes, i.e. the topological tran¬ 
sition happens when m = \/N. This manifests itself in 
Fig.i as the point where two momenta become imagi¬ 
nary. 

Conclusions and outlook: In conformal field theories 
the N~^ term in energy is universal and o nly de pends on 
the central charge of the Virasoro algebraP^IHI. Here we 
find that in the case of topological phase transitions this 
term naturally extends into a scaling function, depend¬ 
ing only on the ratio of the system size to the correlation 


length. Furthermore this scaling function is universal 
for all topologically non-trivial classes of non-interacting 
fermions in one spatial dimension. While the scaling 
function for energy appears to be sensitive to the topo¬ 
logical nature of the transition, this is by no means the 
common situation. For example, the finite-size scaling 
function of the entanglement entropy away from the crit¬ 
ical pointP^, appears to be symmetric across the topo¬ 
logical transitioiP^ (there is still an asymmetric size- 
independent boundary term). 

It is natural to ask whether the scaling behavior 
changes with interactions, especially for models with cen¬ 
tral charge different from c = 1 and c = 1/2 considered 
here. Another direction to explore is relation of the scal¬ 
ing function to the theory of integrable system^^. In 
particular if it may be expressed in terms of solutions 
of Painleve equations, as it happens in e.g. the Ising 
modeP^. As mentioned above, the non-analytic contri¬ 
butions to the finite-size scaling near w = 0 are related 
to the monodromy, i.e. discontinuity across the branch 
cut terminating at w = 0, which happens to be simply 
a second order polynomial in w. An open question is if 
the full /(r(;)-function may be recovered from the mon¬ 
odromy data, specified for all of its branch cuts, through 
the solution of a Riemann-Hilbert probleirP^. 

We are grateful to A. Abanov and I. Gruzberg for valu¬ 
able discussions. The work was supported by NSF grant 
DMR1306734. 
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I. SUPPLEMENTAL MATERIAL TO 
UNIVERSAL FINITE-SIZE SCALING AROUND 
TOPOLOGICAL QUANTUM PHASE 
TRANSITIONS 

a. Different symmetry classes: In this section we in¬ 
troduce model systems for all topologically non-trivial 
symmetry classes in one dimensiorP. As men tioned in 
the main text, the Su-Schrieffer-Heeger modeP^*20| be¬ 
longs in AIII symmetry class, it is a bipartite fermionic 
model with tight-binding Hamiltonian 

N N-1 

'Haiii = '^ticl^jCBJ + ^ t2C^g jCA,j+l + h.c. (13) 
1=1 1=1 

and dispersion relation e{k) = ~ 2 tit 2 cos k. 

For ti t 2 neighboring sites form dimers, where for > 
t 2 all sites are part of a dimer, but for ti < t 2 the two 
sites at the ends of the chain are unpaired. These are 
the two distinct phases with topological index Z = Q or 
Z = 1 respectively. The central charge is c = 1. 

The spinless Kitaev chairP^IlSl falls into symmetry class 
BDI. The Hamiltonian is 

N N-1 

y-BDI = - X! (^^1 + 1 + ^^l^ + l) + 

1=1 1=1 

(14) 

where /r is the chemical potential, t a hopping parameter 
and A the pair creation amplitude. For \fi\ < t the system 
is topological and Majorana edge modes appear, the gap 
size gives m = t — |/x|. The central charge is c = 1/2. 

A generalization to spin-1 fermions with p-wave pair¬ 
ing falls into symmetry class DIlP^l 

N N-1 

'Hdiii = -2Ai ^ c] -b t c] (15) 

N-1 

+ E (^i^tEi.t + ^i^tsVid) 

N-1 

+A 2 Y +^-c- 

0 ’ 7 ^ 0 ’M = 1 

Here /i is the on-site potential, t hopping parameter and 
Ai, A 2 p-wave pairing fields. For \fi\ < t the system is in 
a topological state, where the gap size gives m = t — \yt\. 
Due to Kramer’s degeneracy all Majorana modes appear 
in pairs, thus c = 1. 

Introducing an additional Zeeman splitting of the on¬ 
site potential breaks time reversal symmetry and gives 
a system in D symmetry class. Compared to the Dili 


Hamiltonian we replace the on-site term as follows: 

N N 

no ■ -2m E -2E^ic].T^Ft + M2cJ-^Cyf. 

j=i 

( 16 ) 

The system is topological if t € (mi, /i 2 ), so there are two 
topological critical points with c = 1/2. Focussing on the 
latter the size of the gap yields m = t — yi 2 - 

The CII symmetry class may be represented by a bi¬ 
partite tight-binding model for spin-4 fermions. 

N 

T~^cii — M ^ (4'^) 

a',i=l 

o’,I = l 

t2 

+ 4,N+h1- - + h.c., 

where m is the on-site potential and ti,t 2 are hopping 
parameters. The critical point is at |/i| = ti, with m = 
— M, and the transition is described by a conformal field 
theory with c = 1. 



FIG. 3: Integration contour in the complex k plane for 
A = 10 and w = 1.1. The contour (red) encloses all zeroes 
of the argument of the logarithm (blue) and is to be taken 
counter-clockwise. After a change of variable z = Nk and 
taking the scaling limit A —^ 00 the contour runs along 
the entire real line, above and below. These two parts of 
the contour are now deformed (purple) to run along the 
branch cuts of \/m^ + (green). 

b. Derivation of the scaling function In this sec¬ 
tion we want to provide more detail on the derivation 
of the scaling function ([^. The total energy of the 
filled lower band is the sum of all states, satisfying 
cos{kN + 6{k)) = 0; tani5(A:) = Eq. (7) in 

the main text, substituted in the lower band dispersion 
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relation e“(fc) = +1^ — 2tit2 cos k. This sum may 

be converted into an integral with the argument princi¬ 
ple: 


E{N, m) 


1 

2 


dk 

2Tri 


(/c) In [cos(-/V/c-I-(5(A:))], (18) 


where the ± sign is chosen differently above and below 
the real line so that the remaining exponential term van¬ 
ishes as k ^ ±ioo. Now one subtracts Ne -I- b from 


Eq. (18) to find the scaling function f{w), 


fjw) 

N 


= ± 


dk 


■(/c)i9feln 


1 -I- e 


:f2i{Nk+5{k)) 


■ ( 20 ) 


where the contour runs in the complex fc-plane encircling 
all solutions of the quantization condition, see Fig. 
The prefactor 1/2 is introduced to avoid double count¬ 
ing. The two pieces of the integration contour above and 
below the real axis yield the same contribution, in the 
following we treat these two pieces independently. The 
contributions proportional to and iV° are the bulk 
and boundary terms, respectively. Writing the cosine as 
sum of two exponentials one can factor out these leading 
contributions to get 


N€ + b = J{dk/2n)e-{k)[N + dkS{k)] (19) 

= ±J (dfc/27rf)e-(fc)5fclne=^*(^'=+‘^('=», 


In the scaling limit N ^ oo integration runs along the 
entire real line and can be seen as closed at infinity. The 
approximation to the dispersion relation near the Dirac 
point, e~(fc) = + k^, has a branch cut along the 

imaginary axis starting at fc = ±im. Now one deforms 
the integration contour to run along this branch cut and 
redefines the integration variable as z = ±ikN and intro¬ 
duces w = Nm. As a result, in the scaling limit N ^ oo 
one finds 



( 21 ) 

where 5w{z) = — tanh ^{w/z), which is Eq. (9) of the 
main text. 











